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We consider the role played by the 4/ states in the rare-earth oxyarsenides REOFeAs 
(RE=Ce,Pr,Nd) and the oxyphosphate CeOFeP, using a first-principles technique that combines 
the local density approximation and dynamical mean-field theory (LDA+DMFT). In the Pr and 
Nd compounds, the 4/ states are located well below and above the Fermi level Ef, and essentially 
do not interact with the iron 3d bands located near Ef, resulting in local moment behavior. In 
the Ce compounds, our results reveal a qualitatively different picture, with the 3d-4/ hybridization 
being sufficiently strong to give rise to an observable Kondo screening of the local 4/ moment. Our 
LDA+DMFT electronic structure calculations allow us to estimate the Kondo temperature Tk for 
both CeOFeP and CeOFeAs. For the phosphate, the order of magnitude of our estimate is consistent 
with the experimental observation of Tk '^ lOK. At ambient pressure, Tk is found to be negligibly 
small for CeOFeAs. Under applied hydrostatic pressure, we predict an exponential increase of Tk 
which reaches values comparable to the superconducting Tc —40 K at pressures above lOGPa. We 
conjecture that the competition between the Kondo effect and superconductivity may be at the 
origin of the monotonous decrease of Tc observed in CeOFeAs under pressure. We argue that the 
quantitative aspects of this competition are inconsistent with a weak-coupling BCS description of 
the superconductivity in the oxyarsenides. 
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The discovery of the REOi_a;Fa;FeAs superconduc- 
tors with a transition temperature up to Tc= 55 K 
[l|, I2, y, y, M, la| and the observation of superconductiv- 
ity also in oxygen deficient REFeAsOi_5[7] has triggered 
enormous efforts, both experimentally and theoretically, 
aimed at understanding the electronic properties of rare 
earth iron oxypnictides. The evolution of Tc along the 
series of rare earth elements, and possible correlations 
with structural properties have attracted particular in- 
terest. While LaOi_a;Fa;FeAs has a Tc of 26 K, the crit- 
ical temperature is drastically increased by replacing La 
by other rare earth ions (RE=Ce,Pr,Sm,Nd,Gd), up to 
'^SSK for Sm . Possible correlations between the evolu- 
tion of Tc with changes in the structural parameters due 
to decreasing size of the rare-earth ions along the series 
have been reported in the literature [T], Isj . 

On the theoretical side, most ab initio studies have so 
far concentrated on the La compound. To our knowl- 
edge, only the work of Nekrasov et al. [9] has ad- 
dressed the question of how the electronic structure of 
REOFeAs changes upon rare-earth substitution. On the 
basis of LDA calculations, these authors obtained essen- 
tially identical electronic structures for all the REOFeAs 
compounds considered (with RE=La,Ce,Pr,Nd,Sm,Y). 
In their work [y], the 4/ shells of the rare-earths were 
treated as core states. However, an intriguing question is 
whether the low-energy electronic structure of REOFeAs 
might be modified along the rare earth series, due to the 
hybridization between the Fe 3d band and the locaHzed 
4/ states of the rare earth ions. This possibility has not 



been so far investigated in the literature. Meanwhile, a 
Kondo screening of the Ce local magnetic moment on the 
4/ shell has been observed in the CeFePO compound, an 
homologue of CeFeAsO, with a reported Kondo temper- 
ature Tk = 10 K [10(1 . This impHes the existence of an 
hybridization between the Fe 3d bands in the vicinity 
of Ep and the localized Ce 4/ states. The discovery of 
the Kondo effect in CeFePO raises the possibility that 
the Re4/-Fe3d hybridization might similarly lead to a 
Kondo screening of the RE local moment in some of the 
oxyarsenide compounds. 

In order to investigate the 4/ states and their inter- 
action with other bands, a many-body treatment of the 
strong Coulomb interaction between 4/ electrons is nec- 
essary. Therefore we have employed the combined local 
density approximation and dynamical mean-field theory 
(LDA+DMFT) approach [H IH to study properties of 
4/ states in CeOFeP and REOFeAs (RE=Ce,Pr,Nd). We 
show that in the Pr and Nd compounds the 4f states are 
very localized, leading to unscreen local moments. The 
occupied 4f states are located well below the 3d-Fe states, 
and the unoccupied ones well above. 

We find that the situation is different, however, for 
the Ce compounds CeOFePn (Pn=P and As), for which 
the occupied 4/ band is located at the bottom or just 
below the bottom of the 3d-Fe states. The resulting hy- 
bridization leads to Kondo screening of the 4/ Ce local 
moment at low T. We estimate the Kondo temperatures 
of CeOFeP, and find that the order of magnitude is con- 
sistent with experiments. We predict a rapid increase of 



the (very low) Kondo scale of CeOFeAs as a function of 
pressure, due to the contraction of the Fe-Ce interatomic 
distance and to the corresponding increase of the f-d hy- 
bridization. Moreover, we conjecture that the Kondo ef- 
fect in CeOFeAs may be at the origin of the rather rapid 
suppression of Tc under pressure observed in doped Ce- 
OFeAs [13| (which is in sharp contrast to the behaviour 
under pressure of LaOFeAs [131 at a similar doping level) . 

We have employed the LDA+DMFT approach using 
the recently developed fully-selfconsistent framework de- 
scribed in [1J|. The local self-energy of the 4/ shell is 
computed according to the DMFT prescription and by 
employing the atomic (Hubbard-I) approximation [l5|. 
This approach to local correlations has been shown to 
be appropriate for the localized 4/ shells of rare-earths 
compounds [lj| . We have used the full four- index U ma- 
trix for the local Coulomb interaction. The spin-orbit 
coupHng and splitting of the bare 4/ levels due to the 
crystal field were taken into account. Our calculations 
are performed for the paramagnetic state and at zero 
temperature. 

In order to determine the value of the local Coulomb 
interaction U on the 4/ shell, we have performed con- 
strained LDA calculations. We obtain a value of about 
9.7 eV for both CeOFeAs and CeOFeP. This is substan- 
tially larger than the usual range of U values for pure 
Ce. This result can be explained by the quasi two- 
dimensional environment of the rare-earth sites in the 
case of oxypnictides, with only four nearest- neighbours 
present, which leads to a rather poor screening of the 
local Coulomb interaction. Values of the Slater integrals 
i^^, F^, and F^ , which are known to be weakly dependent 
on crystalline environment, have been taken from the op- 
tical measurements of [l6|. The corresponding values of 
the exchange parameter J are equal to be 0.69, 0.73, 
and 0.77 eV for Ce, Pr, and Nd, respectively. For the 
double counting correction, we have employed the fully- 
localized Hmit expression U{Nf - 1/2) - J{Nf/2 - 1/2), 
where the occupancy Nf in the atomic Hmit and at zero 
temperature is equal to 1, 2, and 3 for Ce, Pr, and Nd, 
respectively. 

All calculations have been carried out at the exper- 
imental lattice parameters a and c [2, l3, lj|- An un- 
usual sensitivity of the low-energy electronic structure of 
iron oxypnictides to the vertical distance between the Fe 
and As planes zas has been pointed out in the literature 
171, ll8|. However, in the present work, we mainly focus 
on high-energy features of the electronic structure, par- 
ticularly on the 4/ bands and their interaction with other 
bands. Those features are not expected to be very sen- 
sitive to small changes in the Fe-As distance, they may 
rather show some sensitivity with respect to the verti- 
cal position of the rare-earth plane zre, which is almost 
constant along the series [s], Il9|, l20||. Thus the experi- 
mental values of the internal parameters zas and zre for 
LaOFeAs [2l| have been used for all oxyarsenides stud- 



ied, while CeOFeP has been computed with the internal 
parameters fixed at their experimental values for PrOFeP 

Hi. 

Because the atomic Hubbard-I approximation to 
strong correlations is not able to capture directly the 
Kondo effect, we have employed the approach of Gun- 
narsson and Schonhammer |23| to calculate the Kondo 
temperature. Within this approach the Kondo temper- 
ature is given (in the Kondo regime) by the following 
expression: 

TK--De "I'^^'^f , (1) 



where e/ is the average position of the occupied 4/ level, 
Po is the density of states at the Fermi level, Nf is the 
/ band degeneracy, Vcf is the average hybridization ma- 
trix element between the 4/ level and the conduction 
band, D is the bandwidth of the occupied part of the 
conduction band. As shown in [23J, the spin-orbit split- 
ting yields a rather small correction to Tk and can be 
neglected, so that we take iV/ = 14. The value of V^f can 
be estimated from the imaginary part of the DMFT hy- 
bridization function at Ep, which reads: 



ImA(£;F) =Im^ 



Ef - (k 



■^PoV^f, 



(2) 



where Ck are the conduction band states. For a degener- 
ate case, the average value of V^r can be extracted from 
Tr Im A (i?i?)/iV/. An accurate evaluation of the parame- 
ters entering in the exponential factor in ^ is necessary 
to obtain any reasonable estimate of Tk- Thus we have 
employed a 9000 k-point mesh in the Brillouin zone in 
our LDA+DMFT calculations in order to evaluated the 
average hybridization Vcf and the 4/ level position e/. 
To obtain an accurate value of the density of states at 
the Fermi level, we use a Full-Potential APW+local or- 
bitals code[2J| with the 4/ states in the core. We have 
checked that all the values used in calculation of Tk are 
converged with respect to the k-point mesh. 

First, we discuss the LDA+DMFT electronic structure 
of the CeOFeP compound in connection with the Kondo 
screening of the Ce local magnetic moments observed ex- 
perimentally in this compound. In the upper panel of 
Fig. [1] the LDA+DMFT spectral function at C/ - J =9 
eV (in red) is superimposed to the band structure ob- 
tained by treating the Ce 4/ states as core (in blue). 
The empty Ce 4/ states (the upper Hubbard band) form 
rather dispersionless bands between 6 and 9 eV above 
Ef hybridized to a certain degree with the Ce 5d band. 
The occupied part (the lower Hubbard band) of the Ce 
4/ states is located at about 1.8-2 eV below Ep- It hy- 
bridizes substantially with the p states of mostly oxy- 
gen character at the top of the As/0 p band, causing a 
downward shift of 0.5 eV of the oxygen 2p band that was 
located at ~ 2.3 eV in the LDA with f in core band struc- 
ture. In the lower panel the corresponding partial density 




band under consideration [25J. Our estimates for e/ and 
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Figure 1: Electronic structure of CeOFeP obtained within 
LDA with the 4/ states treated as core (blue line) and within 
LDA+DMFT (red line). The 4/ states of RE form a rather 
dispersionless upper Hubbard band 6-9 eV above the Fermi 
energy. The lower Hubbard band is more dispersed due to 
hybridization with O 2p and Fe 3d. It is located between -1.8 
and -2 eV/ In the lower panel the density of occupied states 
is displayed, with the total, partial Ce 4/, Fe 3d, pnictogen p 
and O 2p shown by black, red, green, blue, and orange curves, 
respectively. One may notice a rather strong hybridization 
between the O 2p states at the top of p band and the occupied 
4/ states. 



of (occupied) states is depicted. Note that, because the 
Hubbard-I approximation does not capture the Kondo 
peak, the LDA+DMFT (Hubbard-I) spectral function in 
Fig.[T]should be viewed as the band structure for T > Tk- 

From the imaginary part of the Ce 4/ hybridization 
function we have obtained T4/=116 meV, while e/=-2.0 
eV, po = l-85 (eVxformula unit)~^, and D ~2 eV for the 
Fe 2>d bands, resulting in an estimate of Tk — 77K. This 
is a reasonable order of magnitude in comparison to the 



experimental result Tk='\-0 K of Briining et al. |lu|, given 
that our estimate for Tk is obtained for the single impu- 
rity case. Indeed, the corresponding Tk for the Kondo 
lattice is expected to be reduced by at least a factor 
of 2 for the approximately half-filled conduction (Fe 2>d) 



poV^r=25 meV also compare well with the corresponding 
values of -2.4 eV and 19 meV, respectively, obtained for 
CeCu2Si2 by Kang et al. [2a| from a fit to the exper- 
imental PES spectra. As noted in Ref. [10|, CeCu2Si2 
and CeOFeP have about the same Kondo temperature 
and Sommerfeld coefficient 7 of the specific heat. 

We have also performed LDA+DMFT calculations for 
CeOFeP using the "standard" value of 6 eV for the U 
parameter of pure Ce 4/. Those calculations predict the 
lower Hubbard band to be much closer to the Fermi level, 
with e/ ~ — 0.5eV. This would lead to a unreahstically 
high Kondo temperature for CeOFeP, even if one takes 
into account necessary corrections to (U) due to devia- 
tion from the Kondo regime (mixed valence behavior) 
for small values of |e/[. Therefore, the constrained LDA 
value of the local Coulomb interaction U — J=9 eV pro- 
vides a better description of the Kondo behaviour in Ce- 
OFeP as compared to the usual choice of U for pure Ce. 
We used the same value to study the oxyarsenides. 

In Fig. [21 we compare the LDA+DMFT band struc- 
tures of stoichometric CeOFeAs, PrOFeAs and NdOFeAs 
with the band structures obtained for the same com- 
pounds within LDA by treating the rare-earth 4/ states 
as core. Differences between the LDA+DMFT and LDA- 
with-4/-in-core electronic structures are especially pro- 
nounced well above Ep, where the unoccupied 4/ states 
form a set of dispersionless bands in the range of ener- 
gies from 6 to 9 eV, from 4 to 8 eV, and from 3 to 7 eV 
for the Ce, Pr, and Nd compounds, respectively. The 
width of about 3 eV of those upper Hubbard bands is 
due to multiplet splittings associated with the exchange 
and spin-orbit interactions. The occupied 4/ states are 
located in the range of energies from -2.3 to -3 eV in Ce- 
OFeAs, from -4 to -5 eV in PrOFeAs, and below -6 eV 
for NdOFeAs. Similarly to the case of CeOFeP, the occu- 
pied 4/ states hybridize rather strongly with the oxygen 
p states, which are located mainly at the top of the As 
4p/0 2p band. This 4/-2p hybridization leads to some 
modifications at the top of the 0/As p band in compar- 
ison with the /-in-core band structure. In contrast, the 
Fe 3d bands are hardly modified by the interaction with 
the 4/ states, as one may see in the upper panel of Fig. [2] 
where blue and red bands corresponding to Fe (close to 
the Fermi level) essentially coincide. From these results, 
we conclude that the 4/ electrons behave as unscreened 
localized moments in the Pr- and Nd- compounds. 

In contrast, for CeOFeAs, the lower Hubbard band is 
located rather close to the Fe 3d bands (see the corre- 
sponding DOS on Fig. [2]). Thus, in analogy to the case 
of CeOFeP, in CeOFeAs one might expect to observe a 
Kondo screening of the 4/ local moment. We have esti- 
mated the Kondo temperature in CeOFeAs using formula 
(U), and the resulting values are given in the first row of 
Table HI With D ~ 2eV for the Fe3d bands, this leads to 
a value of Tk of order 10~^ K. Thus, Tk in CeOFeAs is 
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Figure 2: Electronic structure of stoichiometric iron oxyarsenides obtained within LDA with the 4/ states treated as core 
(blue line) and within LDA+DMFT (red line). The empty 4/ states of RE form a rather dispersionless upper Hubbard band 
significantly above the Fermi energy. The lower Hubbard band is located just below the bottom of the Fe 3d band, in the 
middle of the 0/As p band and below the p band in CeOFeAs, PrOFeAs, and NdOFeAs, respectively. The lower panel displays 
the corrseponding densities of states, the color coding is the same as in Fig. [T] 



negligible in comparison to the magnetic ordering tem- 
perature for the Ce local moments, which was measured 
to be about 5 -ftT in this compound|4|. This drastic reduc- 
tion of the Kondo scale (by about 5 orders of magnitude 
!) in comparison to Tk for CeOFeP is mainly due to 
a substantially weaker 3d-4/ hybridization in CeOFeAs. 
This decrease of the d-f hybridization is rather expected, 
as the Fe-Ce distance is larger in CeOFeAs than in Ce- 
OFeP. 

In order to investigate the behavior of the Kondo 
scale under applied pressure we have estimated the cor- 
responding change in volume of CeOFeAs. Methods 
based on the atomic sphere approximation can not re- 
liably describe elastic properties of a complex crystal 
structure. One may notice that the structural param- 
eters of CeOFeAs and LaOFeAs are rather similar, while 
the 4/ states of Ce, being still essentially localized, are 
not expected to affect substantially the elastic prop- 
erties of CeOFeAs. Therefore, for CeOFeAs we have 
adopted the theoretical elastic constants computed for 
LaOFeAs [27| within a full-potential approach. With 
those elastic constants the relative values of the lattice 
parameters a/ao(c/co) of CeOFeAs (where a^ and cq are 
the corresponding values at zero pressure) are equal to 
0.985(0.978), 0.970(0.956), and 0.955(0.934) under hy- 
drostatic pressure of 5, 10, and 15 GPa, respectively. 
Anisotropy of the elastic constants leads to a substan- 
tially larger relative contraction for the c- parameter com- 



Table I: The average hybridization Vcf and 4/ level position 
£/, density of states at the Fermi level po, and Kondo tem- 
perature Tk in CeOFeAs as function of applied pressure 
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1.77 
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137.0 



pared with the one for a. 

In Table [J we list, as a function of applied pressure, 
the values of the average hybridization Vcf and 4/-level 
position £/ obtained within LDA+DMFT as well as the 
density of states at the Fermi level. One may see that 
Vcf grows substantially with increasing pressure, while 
the 4/-level position is shifted towards the Fermi level. 
The reduction of the Ce-Fe distance naturally leads to 
an increase in the corresponding hybridization matrix el- 
ements, the contraction in the Ce-0 distance further en- 
hances the strong 2p-Ce 4/ hybridization, leading to 
the upwards shift of e/ . The Fe 3d bandwidth increases 
as well under pressure, leading to a drop in po- However 
the effects due to the sharp rise in Vcf are obviously the 
most important, leading to a reduction of |e/|/(yOoKj) in 
formula |(T|) and to a corresponding (almost exponential) 
increase of the Kondo temperature as a function of pres- 
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Figure 3: Experimental evolution of the superconducting 
transition temperature Ta in CeOo.88Fo.i2FeAs under applied 
pressure[i3|(red dashed curve) and our theoretical estimate 
for the corresponding evolution of the Kondo temperature Tk 
in CeOFeAs (blue solid curve). 



sure (see Tabled. Hence, Tk raises from the negligible 
value of 10~^ K at ambient pressure to a value of order 
100 K under an applied pressure of 15 GPa. 

Within our approach one may only obtain a rather 
crude estimate for Tk- But our qualitative conclu- 
sions should be rather robust: as the pressure in- 
creases the Kondo temperature exhibits an exponential 
growth, from values which are negligible compared to 
the superconducting transition temperature Tc at am- 
bient pressure, up to values of the same order of mag- 
nitude as Tc at pressures above lOGPa. It is interest- 
ing to contrast the predicted enhancement of Tk un- 
der pressure with the behavior of Tc observed experi- 
mentally for CeOo.88Fo.i2FeAs [lal (see Fig. [31). In con- 
trast to LaOa;Fi_j;FeAs at a very similar doping level, 
for which Tc raises with increasing pressure up to 5 
GPa before starting to decrease rather slowly [l3|, in 
CeOo.88Fo.i2FeAs one observes a monotonous and rather 
steep decrease of Tc with increasing pressure. The super- 
conducting transition is actually completely suppressed 
at pressures of about 20 GPa. As one may see in Fig.[3l a 
rapid reduction of Tc sets in at pressures (above 5 GPa), 
at which Tk starts to reach values of order of a few de- 
grees Kelvin. Thus, we may conjecture that the rapid 
suppression of superconductivity observed in the doped 
CeOFeAs under pressure is due to the corresponding 
stabilization of a competing heavy-fermion phase with 
Kondo-screened local moments of the Ce shell. 

The competition between the Kondo effect and s-wave 
superconductivity has been investigated theoretically in 
Ref. [30|, in the framework of the periodic Anderson 
model. This study concluded that with Tk > Tcq/W 
(where Tco is the transition temperature in the absence 
of the Kondo effect and W is the bandwidth of the con- 
duction band), the s-wave superconductivity becomes 
rapidly suppressed, with Tc decreasing exponentially as a 



function of Tk ■ This expresses the competition between 
the Kondo screening energy {r^ Tk) and the condensation 
energy ^ T^q/W of a weak-coupHng BCS superconduc- 
tor. In our case the threshold value of Tk = Tc^/W is of 
order of 0.1 K. However, our estimates of Tk lead us to 
conclude that while the Kondo effect may be responsible 
for a rather rapid suppression of the superconductivity 
in Ce0a;Fi_2;xFeAs under applied pressure, the actual 
superconducting phase is still much more robust than 
what would be predicted by the theory in Ref. [30||- This 
points out to important differences between the nature of 
the superconductivity of the rare-earth iron oxyarsenides 
and that of a weak-coupHng BCS superconductor. 

In conclusion, our study reveals that while for the heav- 
ier rare-earths the electrons in the 4/ shell of rare-earth 
oxyarsenides behave as unscreened local moments, the 
Cerium-based compounds behave in a different manner. 
There, a competition between Kondo screening and su- 
perconductivity takes place under appHed pressure. This 
may be responsible for the suppression of superconduc- 
tivity under applied pressure in the cerium-based com- 
pounds, and suggests that the competing phase at large 
pressure should be a heavy-fermion state with a moderate 
effective-mass enhancement. Orders of magnitude sug- 
gest however that the Kondo effect is less detrimental to 
superconductivity in these compounds than for a weak- 
coupling BCS superconductor, hence providing indirect 
evidence for the unconventional nature of superconduc- 
tivity in the oxyarsenides. This competition between su- 
perconductivity and a heavy-fermion state deserves fur- 
ther theoretical and experimental investigations. 
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